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ORIGINAL ARTICLE 

Genome-wide association analyses suggest NELL1 influences 
adverse metabolic response to HCTZ in African Americans 

JL Del-Aguila 1 , AL Beitelshees 2 , RM Cooper-DeHoff 3 , AB Chapman 4 , JG Gums 3 , K Bailey 5 , Y Gong 3 , ST Turner 5 , 
JA Johnson 3 and E Boerwinkle 1,6 

Hydrochlorothiazide (HCTZ) is one of the most widely prescribed antihypertensive medications. Although it is well known that 
HCTZ is associated with hyperglycemia and hypertriglyceridemia, the mechanisms underlying these adverse effects are not well 
understood. We performed a genome-wide association study and meta-analysis of the change in fasting plasma glucose and 
triglycerides in response to HCTZ from two different clinical trials: the Pharmacogenomic Evaluation of Antihypertensive Responses 
and the Genetic Epidemiology of Responses to Antihypertensive studies. Two single-nucleotide polymorphisms (rs1 2279250 and 
rs4319515 (^ = 0.73)), located at 1 1 pi 5.1 in the NELL1 gene, achieved genome-wide significance for association with change in 
fasting plasma triglycerides in African Americans, whereby each variant allele was associated with a 28mgdl _1 increase in the 
change in triglycerides. NELL! encodes a cytoplasmic protein that contains epidermal growth factor-like repeats and has been 
shown to represses adipogenic differentiation. These findings may represent a novel mechanism underlying HCTZ-induced adverse 
metabolic effects. 
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INTRODUCTION 

Hypertension is major risk factor for stroke, cardiovascular 
and kidney diseases. 1 Hydrochlorothiazide (HCTZ) is a widely 
prescribed thiazide diuretic for initial and combination drug 
treatment of hypertension 2 with over 110 million prescriptions 
annually in the United States, but its ability to produce a variety of 
adverse metabolic effects (AMEs), such as hyperglycemia 3-5 and 
dyslipidemia, 3,6-8 is well known. Despite being recognized for 
decades, the mechanisms of these thiazide-induced AMEs are not 
well understood. The ability of thiazides to increase glucose has 
primarily been attributed to their ability to cause hypokalemia, 8,9 
thereby impairing insulin secretion. However, the associations 
between glucose increase and potassium lowering are mostly 
based on study-level associations of average changes, and we 
previously did not observe a relationship between patient-level 
changes in glucose and potassium, raising questions about this 
mechanism. 10 More recently, additional mechanisms have been 
suggested including their ability to increase visceral and hepatic 
fat accumulation, which could promote insulin resistance. 11 The 
mechanisms underlying thiazide-induced increases in triglycerides 
are even less well understood. It has been suggested that 
thiazides' ability to cause volume depletion produces hemocon- 
centration of lipoproteins. Volume contraction stimulates 
the renin-angiotensin-aldosterone system thereby stimulating 
catecholamine release and subsequent adipose tissue lipolysis. 
However, various studies have found that the degree of hemocon- 
centration is inconsistent. 12,13 Another possible mechanism is 
an incremental reduction of lipoprotein lipase activity due to 



interference in the production, release or action of insulin by 
HCTZ. 13 One way to help identify underlying mechanisms of these 
adverse effects and to distinguish among these competing 
hypotheses is to identify genetic variations that are predictive of 
inter-individual variation in AME after HCTZ treatment. Using a 
genome-wide association study (GWAS) it is possible to identify 
genes or genomic regions associated with a phenotype without 
a priori knowledge of a biologic effect of the gene on the trait 
of interest. In this study, we report the results of the first GWAS 
of HCTZ-induced AMEs in European- and African-American 
hypertensive patients from two different pharmacogenetic 
studies. 



MATERIALS AND METHODS 

Study population 

Phenotype and genotype data were collected from the Pharmacogenomic 
Evaluation of Antihypertensive Responses (PEAR — clinicaltrials.gov identi- 
fier NCT00246519) study and the Genetic Epidemiology of Responses to 
Antihypertensive (GERA— clinicaltrials.gov identifier NCT00005520) study. 
In both studies, participants had mild-to-moderate primary hypertension 
without a history of heart disease or diabetes mellitus. Primary 
hypertension was defined as blood pressure (BP) levels > 140/90 mm Hg 
or current use of prescription antihypertensive medications in the absence 
of a known cause for elevated BP. If at any time during the study protocols 
described below the average diastolic BP (DBP) rose to >110mmHg or 
the systolic BP to > 180 mm Hg, participants were withdrawn from further 
study participation and prescribed effective antihypertensive drug therapy. 
Details about the study designs and exclusion and inclusion criteria have 
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been described previously. 14,15 Briefly, in the PEAR study, individuals of any 
race-ethnicity and gender combination from age 17 to 65 years old with 
mild-to-moderate primary hypertension were recruited. Participants were 
enrolled in Gainesville, Fl, USA; Atlanta, GA, USA; and Rochester, MN, USA. 
All participants were newly diagnosed hypertensives, untreated 
hypertensives or treated hypertensives taking less than three 
antihypertensive drugs. The protocol of the study was as follows: a wash 
out period of approximately 4 weeks was done in order to remove 
the effects of previous BP medication (if any) from the participants. If 
at the end of the wash out period the average seated home DBP 
was>85mmHg, office DBP was>90mmHg and the home and office 
systolic BP was <180mmHg the individuals were enrolled into the 
randomized phase. All biological samples (blood and urine) were collected 
in the fasting state. As soon as the baseline evaluations were completed, 
the individuals were randomized to HCTZ (thiazide diuretic 12.5mg orally 
once daily) or atenolol (^-blocker 50 mg orally once daily) for 3 weeks, with 
dose doubling (25 mg) for those with BP > 1 20/70 mm Hg for an additional 
6 weeks. More than 90% of PEAR participants received the higher 25 mg 
dose of HCTZ. For this AME GWAS analysis, we used only the patients 
randomized to the HCTZ from the PEAR study, referred to hereafter as 
PEAR HCTZ monotherapy. 

In the GERA study, African Americans and European Americans with 
primary hypertension were recruited at Emory University in Atlanta, GA, 
USA and at the Mayo Clinic in Rochester, MN, USA, respectively. The 
protocol for the GERA study was as follows: a wash out period of at least 4 
weeks was done in order to remove the effects of previous BP medication. 
If at the end of the wash out period the average office DBP was 
>90mmHg, qualifying individuals were treated with HCTZ (25 mg orally 
once daily) for 4 weeks. BP was measured in the seated position using a 
mercury sphygmomanometer and blood samples were obtained for 
baseline biochemical measurements. At the end of the 4-week diuretic 
treatment period, BP was measured and blood samples were again 
obtained for biochemical measurements. All blood collections were done 
in the morning after 8h of fasting. 

All patients enrolled in PEAR and GERA provided written informed 
consent, and the institutional review boards of participating study centers 
approved the study protocols. 



Phenotype and genotype data 

All biomedical measurements were made in a central laboratory at the 
Mayo Clinic. In the PEAR study, these methods were implemented on a 
Hitachi 91 1 Chemistry Analyzer (Roche Diagnostics, Indianapolis, IN, USA). 
In the GERA study, plasma glucose concentrations were determined by 
automated spectrophotometric methods implemented on an IL Monarch 
Chemistry system 760 (Instrumentation Laboratories, Lexington, MA, USA). 
Triglyceride concentrations were determined spectrophotometrically using 
Roche reagents on a Cobas Mira analyzer (Roche Diagnostics, Somerville, 
NJ, USA). In both studies, plasma glucose and triglyceride concentrations 
were measured at the end of the wash out period (baseline) and at the end 
of the HCTZ period (final). Glucose and triglyceride responses to HCTZ 
were defined as the difference between the levels at the final and the 
baseline visits. Individuals with response values under or over 3 standard 
errors from the mean response were removed from the analysis. Plasma 
insulin, which was used as a covariate during statistical analyses, was 
measured using the Access Ultrasensitive Insulin immunoassay system 
(Beckman Instruments, Chaska, MN, USA). 

In the PEAR study, individuals were genotyped on the lllumina 
HumanOmnil-Quad (lllumina, San Diego, CA, USA). GERA participants 
from the opposite extremes of the DBP response distribution were 
genotyped using the GeneChip Human Mapping 500k Array (Affymetrics, 
Santa Clara, CA, USA) using standard procedures. As part of routine quality 
control steps, single-nucleotide polymorphisms (SNPs) with minor allele 
frequency <1%, call rates <95%, Hardy-Weinberg equilibrium P-values 
^10 -5 and individuals with >10% missing genotypes were removed 
from the analysis. The software MACH 16 (www.sph.umich.edu/csg/ 
abecasis/MACH/download) was used to impute the approximate 2.5 
million HapMap SNPs using the phase II CEU as the reference panel for 
European Americans and a cosmopolitan sample of CEU and YRI for African 
Americans. Quality control for the imputed results was done using 
standard procedures (RSQ_HAT <0.3 and minor allele frequency ^0.05). 
MACH generated a file with the highest posterior probabilities for each 
imputed SNP, which was used in the analysis. After quality control and 
imputation, there were > 2 million SNPs used for the genotype-phenotype 
association studies in each race-study group. 



Statistical analysis 

Analyses were done by self-reported race within each study and later the 
results were combined across studies within races using meta-analysis 
using a fixed effect model 17 weighted by the inverse of the variance of the 
race and study-specific p-estimates. Linear regression assuming an additive 
model was performed using ProbABEL 18 to evaluate the association 
between each SNP and metabolic responses to HCTZ. The variables sex, 
age, waist circumference, baseline glucose or triglycerides, and baseline 
insulin were used as covariates. The inclusion of the baseline value as 
covariate has both detractors 19,20 and supporters, 21,22 and our goal was to 
identify loci that influence change in glucose and triglycerides 
independent of any baseline effects. 

To avoid the possibility of spurious association as the result of 
population substructure, 23,24 the two first principal component scores 
were used as covariates in the analysis. These principal component 
scores were obtained using the EIGENSTRAT 24 software (genepath.med. 
harvard.edu/~reich/EIGENSTRAT.htm). After meta-analysis, the definition 
of statistical significance was P<5x10 -8 and the same direction of 
effects (that is, the sign of f>) in each study. A statistically suggestive 
P-value was defined as P<1 x 10 -5 and the same direction of effects 
and was used in order to avoid discounting true positive associations. 



RESULTS 

Baseline characteristics of the European-American and African- 
American study participants are shown in Tables 1a and 1b, which 
shows the relative similarity in participant characteristics between 
the two studies. In both PEAR and GERA, there were a higher 
percentage of European-American male participants, while for 
African Americans there were a higher percentage of females. 
Average age and body mass index were nearly identical in PEAR 
and GERA. Baseline triglycerides were higher in European 
Americans from GERA compared with the other groups. Response 
to HCTZ treatment are shown in Table 2, and document that HCTZ 
significantly increased glucose and triglyceride in both race 
groups in both studies. 



Table 1. PEAR and GERA sample characteristics measured at baseline 



{a) Variables 


PEAR: 


GERA: 




European 


European 




Americans 


Americans 


Total 


229 


196 


Gender (% men) 


59.8 


57.1 


Age (years) 


49.9 ± 9.4 


48.5 ± 7.3 


Waist circumference (cm) 


94.8 ±13.6 


103.7 ±14.7 


Insulin (|ilUml -1 ) 


9.4 ±8.7 


8.5 ± 4.9 


Plasma glucose (mgdl -1 ) 


92.9 ±12.2 


92.4 ±12.5 


Plasma triglyceride (mgdl -1 ) 


149.6 ±100.9 


184.3 ±97.55 


Systolic blood pressure (mmHg) 


151.8±13.5 


1 42.6 ± 12.6 


Diastolic blood pressure (mm Hg) 


97.9 ± 5.8 


95 ± 5.7 


{b) Variables 


PEAR: African 


GERA: African 




Americans 


Americans 


Total 


148 


194 


Gender (% men) 


37.84 


48.5 


Age (years) 


47.4 ± 8.9 


47.4 ±6.3 


Waist circumference (cm) 


97.3 ±13.5 


96.8 ±14.4 


Insulin (|ilUml -1 ) 


9.9 ± 9.6 


12.9 ±10.1 


Plasma glucose (mgdl -1 ) 


90.8 ±12.7 


97.4 ±14.4 


Plasma triglyceride (mgdl -1 ) 


105.84 ±93.98 


122.42 ±57.52 


Systolic blood pressure (mmHg) 


151.8 ± 13.14 


150.8 ±14.9 


Diastolic blood pressure 


97.9 ± 5.8 


97.2 ± 5.5 



(mm Hg) 

Abbreviations: GERA, Genetic Epidemiology of Responses to 
Antihypertensive; PEAR, Pharmacogenomic Evaluation of Antihypertensive 
Responses. 
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Table 2. HCTZ treatment response in PEAR and GERA sample 



{a) Variables PEAR: European Americans GERA: European Americans 



Plasma glucose response (mgdl" 1 ) 1.44 ±8.9 (P = 0.01) 3.03 ± 7.6 (P<0.0001) 

Plasma triglyceride response (mgdl" 1 ) 8.21 ±58.7 (P = 0.2) 14.57±59.9 (P<0.0001) 



{b) Variables PEAR: African Americans GERA: African Americans 

Plasma glucose response (mgdl -1 ) 2.55±10.8 (P<0.005) 4.3±9.5 (P<0.0001) 

Plasma triglyceride response (mgdl -1 ) 7.95±40.6 (P = 0.02) 9.56±40.3 (P<0.0001) 

Abbreviations: GERA, Genetic Epidemiology of Responses to Antihypertensive; PEAR, Pharmacogenomic Evaluation of Antihypertensive Responses. 
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Figure 1. Genome-wide association study (GWAS) meta-analysis Manhattan plot showing the association of triglyceride response to 
hydrochlorothiazide (HCTZ) treatment in African Americans in the Pharmacogenomic Evaluation of Antihypertensive Responses (PEAR) and 
the Genetic Epidemiology of Responses to Antihypertensive (GERA) studies. Two single nucleotide polymorphisms (SNPs) on chromosome 1 1 
showed genome-wide significance (P = 5 x 10 ). Eighteen loci showed P-values less than the suggestive threshold (P= 1 x 10~ 5 ). 



GWAS of adverse effect of HCTZ on glucose 
The meta-analysis of glucose change during HCTZ treatment 
revealed no SNPs that achieved the a priori definition of 
genome-wide significance in either European Americans or 
African Americans. The Manhattan plots and QQ plots for these 
analyses are shown in the online Supplementary Figures SI and S3 
(Supplementary Information) However, in European Americans 
48 SNPs from 8 genomic regions achieved our definition of 
suggestive significance with P<1 x 10~ 5 having effects in the 
same direction in both studies. Similarly, in African Americans 
61 SNPs in 12 genomic regions met these criteria. Each of these 
SNPs is summarized in Supplementary Tables S1.1, S1.2 and S3.1, 
S3.2 (Supplementary Information). 



GWAS of adverse effect of HCTZ on triglycerides 
Among African Americans, two SNPs (rsl 2279250 and rs4319515 
(r 2 = 0.73)) on chromosome 1 1 , reached genome-wide significance 
(jff = 28mgdr 1 ,P=6.6x 10" 9 and £ = 27 mgdl" \P=4.05 x 10~ 8 , 
respectively). The Manhattan plot is shown in Figure 1, and 
Figure 2 displays the triglyceride response by genotype in both 
PEAR and GERA. These SNPs are in the NELL1 gene, whose 
encoded protein is involved in adipose cell differentiation. 25 No 
SNPs achieved genome-wide significance for triglyceride response 
among European Americans. 



As was the case with the HCTZ-induced change in glucose, the 
triglyceride response also had many SNPs that met our suggestive 
significance threshold. In European Americans, 25 SNPs in 
10 genomic regions met the threshold of P<1 x 10~ 5 with a 
same direction effect, with 77 SNPs in 27 genomic regions 
meeting this threshold in African Americans. Many of these SNPs 
were genes previously associated with metabolic syndrome, 
diabetes, insulin and other metabolic traits. These SNPs are 
described in detail in Supplementary Tables S2.1, S2.2 and S4.1, 
S4.2 (Supplementary Information). 



DISCUSSION 

We conducted a GWAS and meta-analysis of adverse metabolic 
responses to HCTZ from two studies of hypertensive patients. We 
focused on change in glucose and triglycerides, which are known 
adverse effects of thiazide diuretics and may partially offset the 
antihypertensive benefit of HCTZ. 3 Among African Americans, two 
SNPs in NELL1 achieved genome-wide significance for their 
association with triglyceride response. NELL-1 protein has 
osteoinductive properties 26 and has been found to repress 
adipogenic differentiation in both unipotent preadipocytes and 
in multipotent adipose-derived stromal cells 25 The regulation of 
adipose stores through differentiation is tightly controlled, with 
adverse metabolic consequences of disordered fat storage 
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Figure 2. Means change in triglycerides adjusted for sex, age, waist 
circumference, base line insulin, base line triglycerides, and the two 
principal components depending on rs1 2279250 genotype (CC, CT, 
TT) for Pharmacogenomic Evaluation of Antihypertensive Responses 
(PEAR) and the Genetic Epidemiology of Responses to Antihyper- 
tensive (GERA) studies. Error bars represent s.e.m. 

occurring in the setting of either too much (obesity) or too little fat 
(lipodystrophy). 27 Therefore, we conjecture that HCTZ could be 
modulating adipocyte differentiation through NELL1 leading 
to accumulation of plasma triglycerides in susceptible patients. 
Further studies to elucidate how HCTZ might influence NELL1 
function are necessary. As our GWAS association was obtained 
through a meta-analysis, confirmation of these results will be 
required in an independent sample of hypertensive patients 
treated with HCTZ. However, this effect was observed in two 
independent samples included in the meta-analysis, and the 
magnitude of the effect is clinically relevant. 

In addition to the SNPs achieving genome-wide significance for 
triglycerides, there were numerous other SNPs in both African 
Americans and Europeans Americans that met our suggestive 
threshold of 1 x 10~ 5 with same direction effects. These included 
SNPs in multiple genes that have previously been associated with 
a variety of relevant metabolic traits. For example, in European 
Americans SNP rs1 7560407 in chromosome 5 had the lowest 
P-value (P = 5.41 x 10~ 7 ) for triglyceride response. This SNP is 
close to the myocyte enhancer factor 2C gene, which regulates 
GLUT4 expression in both skeletal muscle and adipose tissue 28 
and has been associated with type 2 diabetes. 29 Another gene 
related to GLUT4 translocation in adipocytes 30 in which we had a 
suggestive signal for triglyceride response is the pleckstrin 
homology-like domain, family B, member 1(PHLDB1) located on 
chromosome 11 (rsl 1216831, P = 4.1 1 x 10~ 6 ). On chromosome 
7, rs4722750 (P = 8.51 x 10~ 6 ) is close to the JAZF zinc-finger 
1 gene, 31,32 which was previously identified in a GWAS of type 2 
diabetes. Similarly, among African Americans we also identified a 
number of SNPs in genes previously associated with metabolic 
traits. For example, genes that have been related to the metabolic 
syndrome include angiopoietin-like 6 (ref. 33) on chromosome 
19 (rs2116940, P = 3.39 x 10~ 6 ), growth differentiation factor 3 
(ref. 34) on chromosome 12 (rs12307997, P = 6.24x10~ 7 ) and 
ubiquitin-like 5 (ref. 35) on chromosome 19 (rs21 16940, 
P = 3.39 x 10~ 6 ). On chromosome 1 and close to rsl 1810574 
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(P = 3.60x 10~ 6 ) is the gene cholinergic receptor muscarinic 3, 
which has been associated with insulin secretion. 36 

For glucose response, no SNP achieved genome-wide signifi- 
cance in either race group, but as with triglycerides, there were 
many reaching our suggestive threshold that have previously 
been associated with diabetes or other metabolic traits. Specifi- 
cally, the region near rs7762018 was associated with type 1 
diabetes in an affected sib-pair analysis, 37 and confirmed later in a 
larger family study 38 and later in GWA studies. 39 Previous meta- 
analysis 40 and candidate gene study 41 have shown a significantly 
increased risk of type 2 diabetes with the polymorphism 
Thr394Thr(G/A) in the peroxisome proliferator-activated receptor 
gamma, coactivator 1 alpha gene in Indians and Asians, which is 
close to rs151 1453 (P = 4.71 x 10~ 6 ) on chromosome 4. The SNP 
rs6859974 is close to the mannosyl (alpha-1,3)-glycoprotein beta- 
1,2-A/-acetylglucosaminyltransferase gene (MGAT1) which is 
related to fatty acid metabolism. 42 MGAT1 was also previously 
associated with obesity in linkage and GWAS. 43 SNPs with 
suggestive associations in African Americans for glucose response 
included SNPs from two genes previously reported in other GWAS of 
T2D, 44 major facilitator superfamily domain containing 9 close to 
rsl 3402330 (P = 5.58 x 10~ 6 ) on chromosome 2 and branched 
chain amino-acid transaminase 1, close to rs7965364 (P = 6x 10~ 6 ) 
on chromosome 12. In another region on chromosome 12, the SNP 
rs7964748 (P=2.45 x 10~ 6 ) is next to the insulin-like growth factor 
1 gene, which is associated to fasting insulin and insulin resistance 
(homeostatic model assessment-insulin resistance) 45 SNP rs225675, 
located on chromosome 6, and is near G protein-coupled receptor 
126 gene, which was associated with body mass index in a previous 
GWAS. 46 

It has been hypothesized that the mechanism by which HCTZ 
leads to hyperglycemia involves increased potassium excretion 
and hypokalemia leading to defects in insulin secretion. 47 
However, this hypothesis has been controversial given that 
many genetic causes of chronic hypokalemia (for example, 
Gitelman's syndrome) are not associated with hyperglycemia or 
increased risk of diabetes 48 Therefore, it is interesting to note that 
none of the genes identified to be suggestively associated with 
change in glucose were obviously related to potassium transport 
or excretion 8 As such, our results are more consistent with a 
mechanism involving insulin resistance, perhaps secondary to 
visceral fat redistribution and hepatic fat accumulation as 
suggested by Eriksson et al. u 

Herein, we report the first GWAS to our knowledge of HCTZ- 
induced AMEs. The strengths of this study include the population- 
based design, the high quality of genotyping and phenotyping, 
the wash out period of all antihypertensives before initiation of 
HCTZ, as well as detailed follow-up of study participants. At the 
same time, we are aware of the limitation posed by small sample 
sizes, especially to detect true effects size. We are also aware of 
possible biases resulting from inclusion of the baseline value as a 
covariate. 19 Although under-powered, it is interesting to note that 
many of our suggestive regions have been identified through 
previous GWAS results as being related to type 2 diabetes, obesity, 
fatty acid metabolism, lipid metabolism, glucose metabolism or 
body mass index or to other a priori identified biologic candidate 
genes (please see http://hyper.ahajournals.org). However, the 
associations we have identified are novel in that they have 
not previously been related to drug-induced adverse metabolic 
side effects. These findings may help provide insight into 
the mechanisms of antihypertensive-induced AMEs. It will be 
important to replicate these findings identified through meta- 
analysis in independent populations of individuals treated with 
thiazide diuretics. 
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